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ABSTRACT 

We report the discovery of fifteen previously unknown Wolf-Rayet (WR) stars found 
as part of an infrared broad-band study of candidate WR stars in the Galaxy. We 
have derived an empirically-based selection algorithm which has selected ^5000 WR 
candidate stars located within the Galactic Plane drawn from the GLIMPSE (mid- 
infrared) and 2MASS (near-infrared) catalogues. Spectroscopic follow-up of 184 of 
these reveals eleven WN and four WC-type WR stars. Early WC subtypes are ab- 
sent from our sample and none show evidence for circumstellar dust emission. Of the 
candidates which are not WR stars, ~120 displayed hydrogen emission line features 
in their spectra. Spectral features suggest that the majority of these are in fact B 
supergiants/hypergiants, ~40 of these are identified Be/B[e] candidates. 

Here, we present the optical spectra for six of the newly-detected WR stars, and 
the near-infrared spectra for the remaining nine of our sample. With a WR yield rate 
of ^7% and a massive star detection rate of ^65%, initial results suggest that this 
method is one of the most successful means for locating evolved, massive stars in the 
Galaxy. 
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1 INTRODUCTION 

Wolf-Rayet (WR) stars - the chemically evolved descendants 
of the most massive stars (i.e. > 25 Mq) - possess intense 
stellar winds and so continually enrich their local environ- 
ments. Although rare in number, high mass-loss rates during 
the WR stage means that they can completely dominate the 
energetics and chemical enrichment of their local interstellar 
medium (ISM). For example, WR stars comprise only 10% 
of the massive stellar content of NGC 3603, but they are re- 
sponsible for 20% of the ionizing phot ons and ^60% of the 
kinet ic energy released into the ISM l|Crowther fc Dessartl 
1 1998T ) . WR stars have also been suggested as possible pre- 
cursors of core-collapse supernovae of types lb and Ic and 
7-ray bursts. 

Our Galaxy provides an excellent laboratory for study- 
ing massive stars, as we can resolve objects on small scales 
and so hope to achieve sample completeness. However, to 
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date ~300 WR star s have been observed in our Galaxy 
(Ivan der ~75 of which have been reported in 

the last five years. With ~1000-2500 WR stars expected 
to be located with in the Milky Way l|Shara et all Il999l ; 
van der HuchtlliuOll h it is highly likely that a large num- 
ber of WR stars are still waiting to be uncovered. 

Previous WR surveys have utilised optical narrow- 
band interference filters tuned to the strong, characteris- 
tic A4686A Heii/C ill emission feature. These have proved 
to be extremely successful, with the la st large-scale opt i- 
cal survey reporting 31 new WR stars l|Shara et alJll999h . 
However, the Galactic Plane contains large amounts of gas 
and dust, such that high extinction means that many WR 
stars may remain "hidden" from view. Even if optical sur- 
veys were to go as deep as V=25 mag, it is estimated that 
one-third of the Galaxy's WR population would still remain 
undetected (|Shara et al.l fl999h As a result one must turn to 
longer, infrared (IR) wavelengths in an attempt to reveal this 
"missing" population. IR surveys have generally followed a 
narrow-band imaging approach in order to detect WR stars. 
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WR features are much weaker at near-IR wavelengths, but 
for locating WR stars in ste llar clusters this met hod has 
proved to be very successful (jCrowther et al.ll2006h . Unfor- 
tunately, for large scale surveys, especially those concen- 
trated towards the Galactic Centre, the observed success 
rate is very low. iHomeier et al] l|2003h confirmed only 4 new 
WR stars found using IR narrow-band images. With a lack 
of suitable filters for large-scale survey instruments (e.g., 
WFCAM on UKIRT), such surveys require large amounts 
of telescope time and produce large amounts of data. With 
many known WR stars found to be isolated in the field, it 
would appear that a significant number of WR stars may 

still be waiting to be found. 

IWilliams fc Antonopouloul ljl98l| ) and IVan Dvk et al.1 
(2006a|) have showed that it is possible to distinguish WR 
stars via their broad-band near-IR colours. A combination 
of continuum and line emission properties of WR stars dis- 
tinguishes the m from other stellar po pulations at IR wave- 
lengths, but as I Van Dvk et all (|2006al ) have found, they are 
not so distinct that a criterion devised for candidate selec- 
tion would lead to an acceptable (in terms of telescope time) 
predicted success rate. After confirming these suspicions in 
a pilot spectroscopic observing run, they incorporated MSX 
8.0/im data to their colour selection in order to broaden the 
colour baseline. This resulted in the discovery of one new 
WR star suggesting that a combination of near- and mid-IR 
data can potentially provide a better discriminant for WR 
stars. 

Taking this scheme one step further, we have devel- 
oped broad-band, IR colour selection criteria which ex- 
ploits 2MASS (near-IR) and the mid-IR Spitzer GLIMPSE 
(Galactic Legacy Infr ared Mid-Plane Survey Extraordinaire, 
iBeniamin et alj|2003h surveyfl GLIMPSE imaged the inner 
regions of the Milky Way within ±1° of the Galactic Plane at 
3.6, 4.5, 5.8 and 8.0 /im with a sensitivity and resolution that 
excellently matches that of 2MASS ( for a more detailed de - 
scription of the GLIMPSE project see lBeniamin et aDl2003h . 
Covering only 220 square degrees of the Galactic Plane, the 
merging of GLIMPSE and 2MASS colours excludes the pos- 
sibility of an all plane survey, but well complements previous 
WR searches, allowing progress potentially to be made in the 
search for "hidden" WR stars in our Galaxy. 

This paper presents the discovery of 15 new WR stars 
confirmed via either optical or near-IR spectroscopic follow- 
up obs^rvatons^t cornplements a similar study conducted 
bv lMarston et al] (|2006l . in prep) who report the discovery 
of six new WR stars. This paper is structured as follows: 
Sections 2 and 3 briefly summarise the empirical colour se- 
lection algorithm used for target selection and information 
regarding spectroscopic follow-up observations. The newly 
discovered WR stars are presented in Section 4, where we 
consider spectral classification and provide reddening and 
distance estimates for each star. Although the primary goal 
of this survey was the detection of WR stars, we have identi- 
fied a large number of emission-line objects. We briefly dis- 
cuss the non-WR detections in Section [6] Finally, we sum- 
marise our results and draw conclusions in Section [7] 



1 GLIMPSE enhanced products are available at 
http: / / data.spitzer.caltech.edu/popular /glimpse/ 



2 CANDIDATE SELECTION 

In the following section we present a brief outline of the 
selection criteria used to generate an initial candidate list. 

At IR wavelengths WR stars exhi bit free-free emis- 
sion c haracteristic of a dense stellar wind (jWright fc Barlowl 
1975). This emission is superimposed upon the stellar con- 
tinuum such that WR stars exhibit a flatter spectral energy 
distribution relative to normal early-type stars. In a future 
paper we will describe in detail the mechanisms which drive 
the population of WR stars and subclasses to exhibit their 
near- to mid-IR colours, f ollowin g the phenomenology de- 
scribed bv lVan Dvk et al] (|2006bh . 

We have considered the positions of the observed 
colours of known Galactic WR stars in all possible mid- and 
near IR colour-colour combinations (i.e., GLIMPSE only, 
2MASS only and GLIMPSE + 2MASS) to investigate if it 
is possible to distinguish WR stars via their broad-band IR 
colours. We find that the known WR stars are best distin- 
guished in two possible colour-colour combinations. 

(i) Mid-IR colours: Figure QJa) shows the mid-IR 
colours for stars contained in a 1° x 1° degree of the Galac- 
tic Plane centred on Z=312° (6±1°). Overplotted are the ob- 
served mid-IR colours of the known WR stars, which clearly 
separate from the general stellar locus. WR stars appear 
to exhibit a mid-IR excess, with colours [3.6] - [8.0] > 0.5 
and [3.6] - [4.5] > 0.1. WN and WC stars appear to ex- 
hibit similar colours, such that mid-IR colours alone do not 
differentiate between WR subtypes. 

(ii) Mid + near-IR colours: Galactic WR stars also 
appear to be distinguished in a combination of near and 
mid-IR colours. The 2005 "reliable" GLIMPSE catalogue 
was cross-matched with the 2MASS all-sky point source cat- 
alogue using a search radius of 1" to create the colour se- 
quence shown in Figure [Tib). As expected, the majority of 
the objects follow the reddening vector, while the WR stars 
clearly exhibit an 8.0/im excess. The WR stars appear to 
follow a separate, well-defined path relative to the redden- 
ing vector. Again WN and WC stars exhibit similar colours, 
however, WC stars with unresolved circumstellar dust shells 
have a tendency to display the largest 8.0 fim excess. 

Using the positions of the known WR stars we have 
outlined two regions of colour space (as indicated in Fig[JJ 
which define our WR candidates. Of the ~2xl0 7 objects 
common to both to the GLIMPSE and 2MASS catalogues, 
~100 000 were initially selected as candidate WR stars. 

To reduce this number of candidates we have im- 
posed additional selection criteria and restricted ourselves 
to only the most reliable sources, i.e., objects that are un- 
blended, unconfused and of the highest photometric qual- 
ity. This does, of course introduce a sample bias against 
highly crowded regions, i.e., stellar clusters. In addition to 
these "quality flags", we have applied a 2MASS only (J-H 
vs H-Ks) colour selection c riterion using t h e WR near-IR 
colour relation reported by IVan Dvk et al.l l|2006al ) . Addi- 
tionally, ass uming that optica l WR surveys are complete to 
B=14mag l|Shara et al.l 11999). the final selection criterion 
was to place a B>14 magnitude limit for all sources which 
have optical counterparts within a 3" radius listed in the 
USNO-A2.0 catalogue. 

Including all of the above restrictions the final number 
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Figure 1. IR Colour-Colour diagrams for a l°x 1° degree of the Galactic Plane centred on (=312° (fe±l°). Previously reported 
Galactic WR stars (open symbols) clearly separate from the general stellar locus (points), and have been grouped into WN (squares), 
WC (diamonds) and WC9d (circles) stars. The dashed lines indicate the limits used over all GLIMPSE fields to select targets for 
spectroscopic follow-up. For a few cases, we find that their observed IR colours do not fall within the limits of our selection criteria. 
However, all of these are located in complex environments which would suggest unreliable photometry for these stars. 



of WR candidates was reduced to ~5700. Candidates were 
checked against known objects listed in the CDS SIMBAD 
database (http://simbad.u-strasbg.fr/ ), within a 10" radius. 
Objects which already had spectroscopic classifications were 
excluded from our sample (~1%). An additional ~5% were 
listed as possible emission line objects or variable stars. The 
rest had no identification in SIMBAD. 

In order to observe a representative sample for each 
longitude bin, targets were grouped by magnitude, and a 
subsample was selected at random from each magnitude bin. 
Objects that appeared to be extended in 2MASS or Digitized 
Sky Survey (DSS) images were excluded. 



3 SPECTROSCOPIC OBSERVATIONS 

The discovery of the fifteen new WR stars presented here 
is the result of four spectroscopic observing campaigns. 
Optically bright WR candidates in the Galactic latitude 
range of Z=284°-350° were observed from the Cerro Tololo 
Inter-American Observatory (CTIO) on the nights of 2006 
March 19-24. Near-IR spectroscopy of candidates within the 
GLIMPSE range Z=0-65° was obtained from the Apache 
Point Observatory (APO) 2006 June 10-11 and the Infrared 
Telescope Facility (IRTF), Mauna Kea on 2006 August 8-9 
and 11 and September 4-5. For classification purposes we 
have also observed a sample of known WR stars and transi- 
tion stars. An observational log is presented in Table [T] 



3.1 Optical Spectroscopy 

Long-slit spectroscopic observations of 97 of our southern 
hemisphere candidates were obtained at the CTIO 4m- 
Blanco Telescope with the RCSpec spectrograph. Seeing 
conditions were ~0.5"-1.0", and all targets were observed at 
an airmass <1.6. Exposure times ranged between 300 s for 
our brighter targets (B~15mag), to 1800s for the fainter 



objects (B~19mag). The data were obtained using the 
KPGL3 grating with a 1" slit, resulting in a dispersion of 
~1.2 A/pixel and resolution of 3.8 A, as measured from com- 
parison arc lines. The wavelength range of the final spectra 
was ~4000-7500A. 

The data were cleaned of cosmic ray hits, de-biased and 
divided by a normalised flat field. Spectra were extracted 
using standard reduction routines within iraf. Bad rows 
and pixels were masked assuming that they corresponded to 
~10a outliers on flat-field images; any additional flaws were 
identified using the spectrum of the photometric standard 
star LTT 62648. Wavelength calibration arcs were taken at 
regular positions across the sky to allow for instrument flex- 
ure. Atmospheric absorption features have not been removed 
from the final spectra. 

The majority of targets showed evidence of high extinc- 
tion as the signal-to-noise ratio {S/N) of the spectra rapidly 
deteriorated at shorter wavelengths. The blue (~4500A) 
continuum S/N ranged from ~l-30, versus ~30-100 in the 
red (^6800A). Nevertheless, blue WR features were typi- 
cally detected at the 25<r level. 



3.2 Near-IR Spectroscopy 

Near-IR spectroscopic follow-up observations of 87 can- 
didate WR stars were acquired using the ARC 3.5- 
m wit h CorMASS (Corn ell Massachusetts Slit Spectro- 
graph, IWilson et al.ll200lh and the IRTF 3.0-m with SpeX 
jRavner et al.ll2003l ). 

CorMASS and SpeX are both near-IR cross-dispersed 
spectrographs. CorMASS is a low resolution detector and 
has a fixed slit width of 2'.'0. SpeX is a medium resolution 
and was used in the short- wavelength cross-dispersed (SXD) 
mode with a 0'.'8 and 0'.'3 slit, resulting in R=750 and 2000, 
respectively. All near-IR observing modes resulted in a spec- 
tral coverage of 0.8-2.4/im. 
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Table 1. Observing Log 



ID 


Date 


Observations 


Slit 


R 


Wavelength 


Candidates 


WR 








(") 




Coverage 


Observed 


Discoveries 


(1) 


2006 March 19-24 


4.0m-Blanco/RCspec 


1.0 


900 


3900-7500A 


97 


6 


(2) 


2006 June 10-11 


3.5m- ARC/CorM ASS 


2.0 


300 


0.8-2.4jUm 


26 


2 


(3) 


2006 August 8-9, 11 


3.0m-IRTF/SpeX 


0.8 


750 


0.8-2.4/im 


41 


6 


(4) 


2006 September 4-5 


3.0m-IRTF/SpeX 


0.3 


2000 


0.8-2.4^m 


20 


1 



In order to maximise the number of candidates ob- 
served, a limiting magnitude of K^10.5 and K^ll.Omag 
was imposed for Spex and CorMASS spectroscopy, respec- 
tively. 

The Spex and CorMASS data were reduced using the 
Interactive Data Language re duction software too ls Spex- 
tool v3.3 and Cormasstool v3.4 (|Cushing et al.ll2004r ). As the 
data were acquired using the standard "nod" technique, in- 
dividual 2D frames were divided through by a normalised 
flat field before each pair was subtracted to produce a 
background-subtracted image; any residual background was 
subsequently removed during the extraction process. Spec- 
tra were extracted using an appropriate aperture radius, 
which for the majority of sources was 

Atmospheric features were removed from the final spec- 
tra by obser v ing a selection of A0V standard stars following 
IVacca et all (|2003l l. The near-IR spectra of the nine of the 
bona-fide WR stars are shown in Figure IA3I 



4 RESULTS: 15 NEW WR STARS 

Of the 184 candidates spectroscopically observed we have 
confirmed the discovery of 15 previously unknown WR stars. 
For brevity we have assigned a HDM nomenclature to these 
new discoveries. 2MASS designations and magnitudes of 
these new discoveries are presented in Table [2] Kg-band 
finding charts can be found in Appendix iBl 

4.1 Spectral Classification 

Classification of WR stars via their optical spectra is rela- 
tively straightforward as schemes using optical line diagnos- 
tics a re well defined (|Smith et al.lll99q : ICrowther fc Srnithl 
1 19961 ). In the near-IR, classification is somewhat more 
difficult as there are no line dia gnostics which u n iquely 
distinguish all WR subtypes (IFiger fc McLean! 1 19971 ; 
ICrowther et al.|[2006h . Optical and near-IR classification will 
therefore be discussed separately. 

4-1.1 Optical Classification 

The reduced optical spectra (normalised to the continuum) 
are presented in Figure |A"T1 grouped into WN and WC sub- 
types. Prominent WR emission features are identified. In 
order to classify the newly discovered WR stars, we have fit 
Gaussian line profiles to the prominent WR features. Mea- 
sured equivalent widths and FWHM of the prominent opti- 
cal WR features are presented in Table 

As Fig [2] shows, HDM 5 and 6 both display strong 
carbon features associated with late-type WC subtypes. 




Figure 2. Optical CTIO/Rcspec spectra of the prominent yellow 
WC features (left) and blue WN (right) features. 

Following the q uanti tative WC classificatio n scheme of 
ICrowther et al.l (|l998T ). refined from that of ISmith et al.l 
l |l99d ). we assign WC7 and WC9 subtypes to HDM 5 and 
HDM 6, respectively. 

Late-type WC stars have the potential of being sur- 
rounded by circumstellar dust shells which contribute ther- 
mal emission to the stellar continuum at mid-IR wave- 
lengths, depending on dust temperature and density. Using 
optical spectroscopy it is not possible to directly confirm 
the presence of dust. However, both HDM 5 and 6 display 
near-IR colours consistent with non-dusty WC stars. 

The remaining four WR spectra display He and N emis- 
sion features indicative of WN stars. F or this analysis we 
have adopted the classification scheme o f lSmith et ai](|l996T l 
and the observed He I A5876 / Hell A5411 line ratios indi- 
cate a WN7 subtype for HDM 3, and WN6 for both HDM 2 
and 4. 

The spectrum of HDM 1 exhibits distinct differences to 
the other WN stars in our sample. The spectral features are 
notably narrower (FWHM~7A) and lower excitation fea- 
tures are present. We compare the blue (AA4600-5000) and 
yellow (AA4600-5000A) spectrum (normalised to the con- 
tinuum) of HDM 1 with that of Galactic star NS4/WR105 
(WN9) and LMC star Sk-66 40 (WN10) in Figured The 
"blue WR bump" (He n, N ii-iii and C ill) is clearly present, 
as are yellow Nni, Cm and Si ill features. The spectral ap- 
pearance of HDM 1 is remarkably similar to that of the com- 
parison stars, and from that standpoint we classify HDM 1 
as WN9/10. Spectral types WN9-10 are distinguished us- 
ing the relative strengths of the blue N ii-i v emission lines 
jSmith et al.lll996l : ICrowther fc Smithlll997r i. Poor blue S/N 
prevents a classification based on this method for our object, 
so we have used th e HellA4686 vs. the Hel A 5876 lines as a 
subtype indicator (|Crowther fc Smith! 1997). This favours 
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Table 2. Data for the newly discovered Wolf-Rayet stars found by our survey. JHK B magnitudes are taken from the 2MASS catalogue 
and mid-IR magnitudes correspond to values listed in the 2005 reliable GLIMPSE catalogue. Blue and red magnitudes are those listed 
in the 2MASS catalogue and are taken from the USNO-A2.0 Catalogue. 



Star 2MASS I b B R J H K s [3.6] [4.5] [5.4] [8.0] Obs 

Designation ID 



HDM1 


11423766- 


-6241193 


295.13 


-0.86 


18.5 


16.6 


11.50 


10.54 


9.81 


9.07 


8.66 


8.42 


7.60 


(1) 


HDM2 


12065647- 


-6238304 


297.85 


-0.21 


17.8 


15.3 


11.55 


10.82 


10.22 


9.56 


9.19 


8.99 


8.59 


(1) 


HDM3 


12133878- 


-6308580 


298.68 


-0.59 


17.7 


14.5 


10.56 


9.64 


9.03 


8.25 


7.99 


7.74 


7.35 


(1) 


HDM4 


12461614- 


-6257234 


302.34 


-0.09 


19.6 


16.5 


12.18 


11.22 


10.56 


9.81 


9.43 


9.21 


8.27 


(1) 


HDM5 


13101207- 


-6239065 


305.08 


+0.14 


18.1 


15.7 


11.06 


10.09 


9.32 


8.59 


8.32 


8.13 


7.79 


(1) 


HDM6 


16113927- 


-5205458 


331.03 


-0.50 


18.3 


15.3 


10.16 


9.25 


8.46 


7.86 


7.48 


7.36 


6.89 


(1) 


HDM7 


18094505- 


-2017103 


10.23 


-0.41 






12.58 


11.06 


10.14 


9.10 


8.74 


8.49 


8.07 


(3) 


HDM 8 


18131420- 


-1753434 


12.72 


+0.02 


19.2 


15.2 


9.62 


8.60 


9.94 


7.32 


6.87 


6.71 


6.34 


(3) 


HDM9 


18192219- 


-1603123 


15.04 


-0.39 


16.5 


14.3 


9.09 


8.30 


7.76 


7.13 


6.80 


6.65 


6.22 


(3) 


HDM 10 


18250024- 


-1033236 


20.53 


+0.98 






12.21 


11.24 


10.61 


9.74 


9.33 


9.03 


8.67 


(•1) 


HDM 11 


18255310- 


-1328324 


18.05 


-0.57 






10.32 


9.52 


8.96 


8.27 


7.95 


7.79 


7.35 


(3) 


HDM 12 


18400865- 


-0329311 


28.54 


+0.92 






11.77 


10.86 


10.30 


9.54 


9.28 


8.99 


8.64 


(2) 


HDM 13 


18411070- 


-0451270 


27.44 


+0.06 






11.05 


10.17 


9.45 


8.80 


8.32 


8.20 


7.74 


(3) 


HDM 14 


19081797- 


-0829105 


42.40 


+0.15 






10.59 


9.47 


8.71 


8.09 


7.60 


7.38 


6.97 


(2) 


HDM 15 


19334401- 


-1922475 


54.91 


-0.16 


16.3 


13.7 


10.20 


9.61 


9.07 


8.13 


7.75 


7.48 


7.01 


(3) 



Table 3. Emission line equivalent widths (Wa) and FW HM ( both i n A) of prominent op t ical W R emission features. WR subtypes have 
been assigned using the classification schemes of lSmith et al. | jl996r i and lCrowther et al. 



Star Niii/Ciii Hen Hen Cm Civ Hen Ho Spectral 

(4640/50) (4686) (5411) (5696) (5808) (5876) (6560) Type 

FWHM W A FWHM W A FWHM W A FWHM W A FWHM W A FWHM W A FWHM W A 



HDM 1 


9 


10 


9 


10 






7: 


2: 






7 


30 


9 


85 


WN9-10h 


HDM 2 


30 


47 


22 


132 


25 


25 






26 


25 


22 


13 


28 


56 


WN6o 


HDM 3 


20 


17 


14 


33 


11 


5 






weak 


weak 


15 


7 


24 


57 


WN7h 


HDM 4 


29 


41 


22 


118 


21 


30 






27 


28 


21 


10 


26 


65 


WN6o 


HDM 5 


25 


343 


25 


119 


25 


30 


37 


83 


42 


212 


31 


289 






WC7 


HDM 6 


19 


107 


21 


42 


12 


10 


21 


68 


24 


292 


22 


67 






WC9 



lined stars based on the FWHM of the He II spectral features. 
A star is considered to be broad lined if FWHM (Hen 4686 ) 
>30A and/or FWHM (Hen 5411) >4uA (|Smith et al.lll996l ). 
All of the WN stars for which optical spectroscopy is avail- 
able are consistent with weak/narrow lined WN stars. 

WR stars are normally hydrogen deficient objects, so 
the presence of hydrogen in WN stars has an important role 
from an evolutionary standpoint. A detailed abundance de- 
termination would, of course, require a quantitative spectral 
analysis, but it is possible to assess the presence of hydrogen 
using Pickering and Balmer lines at optical wa velengths. 

Using blue Hen / Hen + H/3 line ratios, ISmith et al.l 

l|l996l ) denote the presence of hydrogen with an additional o, 
(h), h, classification to signify no, weak or strong hydrogen. 
For our sample, poor S/N data prevents a reliable analy- 
sis of He II A4541 and so we have qualitatively assessed our 
spectra using Hell + H/3 (A4860), Hell+ Ha (A6560) and 
Hen (A5411). Only for HDM 1 and 3 are the hydrogenic Hell 
lines strong relative to Hell A5411, suggesting a WN9-10h 
and WN7h classification for these objects. 



4-1-2 Near IR Classification 

Recently, with reference t o optical line di a gnosti cs (i.e., he- 
lium and carbon lines), ICrowther et al.l (|2006l ) presented 



+ 
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Figure 3. Comparison of HDM 1, NS4/WR105 (WN9) and LMC 
star Sk-66 40 (WN10). The spectra of NS4/W R1 05 and SK-66 40 
are ta ken fromlBohannan & Crowthcr (1999) and Crowthc r et al.l 
l|l995l l. 



a WN9 subtype for HDM1, but an accurate classification 
requires further observations. 

WN stars can be further classified as weak or broad 
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Figure 4. Spectral comparison of the near-IR C in / C IV mor- 
phology of WR 121 and WR 135 with HDM 13 and 15. All spectra 
were acquired with the IRTF and SpeX. 
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Figure 5. A spectral comparison of the near-IR Nm / N V mor- 
phology of the newly discovered WN7 (left) and WN6 (right) stars 
with that of WR158 (WN7) and WR136 (WN6b). 



a near-IR classification scheme which showed that these 
line diagnostics serve as reliable subtype discriminators. For 
stars earlier than WN6 or WC7, classification requires ad- 
ditional nitrogen and carbon diagn ostics. Here, we have a 
adopted the classification scheme of ICrowther et all (l2006t ) 
as well as making spectral comparisons with previously 
known WR stars. Line strengths of the prominent near-IR 
WR features are summarised in Table [3] 

HDM 13 and 15 both display strong near-IR carbon fea- 
tures and from z and K-band Civ/Cm line ratios we infer 
WC9 and WC8 subtypes for HDM 13 and 15. These classi- 
fications are confirmed in Figure [4] where we compare the 
C IV 1.19 and C in 1.20 /im features of the newly discovered 
WR spectra with those of well- classified WC stars. 

The presence of circumstellar dust is easily inferred from 
near-IR spectra as the WR K-band spectral features are di- 
luted by dust emission. Of the newly presented WC spectra, 
none show evidence of dust emission. 

For HDM 7-12 and HDM 14 spectral features are con- 
sistent with nitrogen-rich WR stars. Hen 1.012/im / He I 
1.083/im and Hen 2.189/im / Br7 line ratios suggest a WN7 
subtype for HDM 8-10 and 12, whereas we infer a WN4-6 
subtype for HDM 7, 11 and 14. 

For WN4-6 stars, He line diagnostics do not provide a 
unique discriminator and the spectral type must be further 
refined using the 2.11 /im Nm / Nv morphology. For the 
stars presented here Nm is strong relative to Nv, and all 
have N line ratios consistent with a WN6 subtype (see Fig[5]). 

In the near-IR, WN stars are classified as broad- 
lined stars if the FWHM (He n 2.1885^m) > 130 A 
ICrowther et all l2006t ). Using this criterion, HDM 8 and 
HDM 10 both qualify as being a broad-lined WN, such that 
we assign WN7b to these stars. 

The presence of hydrogen can be assessed from the near- 
IR spectra using the Hen 2.1889 /im / Br7 line ratios. Ob- 
served line ratios for the WN6 stars mimic those of the H - 
deficient WN6 stars in Westerlund 1 l|Crowther et a.l.ll2006h . 
such that we strongly suspect that our WN6 stars are H- 
deficient. The same conclusions can be drawn for the narrow 
lined WN7 stars. 



4.2 Extinction and Distance 

WR photometry can be used to estimate the extinction 
and distance towards our stars, providing we have a cal- 
ibration of intrinsic colours and absolute magnitudes. In 
principle, both visual and IR wavelengths can be used to 
derive reddenings for our sample, but with extinction at 
JHK wavelengths expected to be relatively insensitive to 
the line-of-sight we have chosen to assess the reddening 
and distances of our star s using their 2MASS JHKs colours 
jlndebetouw et alJliooBT) . 

Interstellar K-band absorption has been estimated by 
comparing the observed and subtype-dependent intrinsic J- 
Ks and H-Ks colours. Heliocentric distances have been de- 
rived using the mean value of Ag g . Th e A[ A j/Ak s ratios 
are taken from llndebetouw et ah (2005) and intrinsic WR 
photo metric properties are those listed in ICrowther et alj 
(2006). Intrinsic WR colours correspond to ESO/SOFI 
JHKs filters whereas absolute magnitudes are based on 
2MASS photometry. Since SOFI and 2MASS have compara- 
ble filter passbands, there should be negligible differences be- 
tween SOFI and 2MASS intrinsic WR colours. Galactrocen- 
tric radii (Rg) have been calculated assuming R Q = 8.5 kpc. 

Derived K-band extinction and distance estimates for 
our WR sample are summarised in Table [S] In general we 
find that Ak s ~0.8-l mag, except for HDM 7 where Ajf s is 
significantly larger. This higher extinction is supported by 
DSS images which shows no evidence of an optical coun- 
terpart. All of the other WR stars in our sample appear 
to have optical counterparts, even though for some optical 
magnitudes are not listed in the USNO-A2.0 catalogue. 

For the majority of our stars, A J t T K and A^J _if are in 
good agreement (to within ~0.1 mag). WR stars are known 
to exhibit a wide range of intrinsic colours, and it is this 
scatter that leads to large uncertainties in our derived ex- 
tinctions. The colour -subtype scatter of the stars which 
ICrowther et alj l|2006l ) used to estimate intrinsic WR colours 
lead to 10% uncertainties in derived extinctions. 

For HDM 5, A J R r K and h^f K differ by 0.4 mag, sug- 
gesting that our adopted WC7 intrinsic colours are un- 
suitable. If this star has intrinsic IR colours more consis- 
tent with those observed for WR68 (i.e., [J — K] o =0.24, 
[H - K] o =0.32, van der Hucht 2006), k J £ K _ and°Ag~ y 
would be consistent to within 0.2 mag, and A« would in- 
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Table 4. Emission line equivalent widths (W\, in A) of near-IR emiss i on lin e diagnostics used to classify the new WR stars. WR subtypes 
have been assigned using the classification scheme of lCrowther et ah (2006).* indicates that it was not possible to de-blend the spectral 
feature and values quoted refer to the strength of the whole feature. 



Star 


He II 


He I 


Nv 


Nni/Heii 


Hell+Br7 


Hen 


Spectral 




(1.012/xm) 


(1.083/im) 


(2.110(um) 


(2.115(um) 


(2.165/Ltm) 


(2.189/xm) 


Type 


HDM7 


221 


187 


40 


90 


7 


31 


WN6o 


HDM8 


132 


338 


10 


45 


59 


62 


WN7b 


HDM9 


54 


244 


52 


32 


3 


38 


WN7o 


HDM 10 


230* 


720* 


9 


65 


87 


77 


WN7b 


HDM11 


194 


207 


42 


84 


6 


30 


WN6o 


HDM 12 


69 


350 


58 


28 


5 


49 


WN7o 


HDM 14 


262 


262 


55 


102 


9 


34 


WN6o 




Cm 


Cn 


Civ 


Cm 


Civ 


Cm 






(0.971/mi) 


(0.990/im) 


(1.191/un) 


(1.198/nm) 


(2.076/im) 


(2.110/Ltm) 




HDM 13 


95 


70 


11 


73 


9 


78 


WC9 


HDM 15 


363 


19 


51 


12 


270* 


95 


WC8 



Table 5. Estimates of the extinction (Ak) and distance (D) de- 
rived using 2MASS near-IR col ours. Intrinsic colours a nd absolute 
Ks magnitudes are taken from lCrowther et al,l 1 120061) . Galactro- 
centric radii (Rg) have been calculated assuming Rq = 8.5 kpc. 



HDM 


SpType 


Kg 


A J-K 




A K 




D 


Rg 


1 


WN9-10h 


9.81 


1.12 


1.04 


1.08 


-5.92 


8.5 


9.1 


2 


WN6o 


10.22 


0.79 


0.76 


0.78 


-4.41 


5.9 


7.8 


3 


WN7h 


9.03 


0.91 


0.93 


0.92 


-5.92 


6.4 


7.8 


4 


WN6o 


10.56 


0.91 


0.97 


0.94 


-4.41 


6.4 


7,1 


5 


WC7 


9.32 


0.36 


0.75 


0.55 


-4.59 


4.7 


7.0 


6 


WC9 


8.46 


0.96 


0.98 


0.97 


-6.30 


5.7 


4.5 


7 


WN6o 


10.14 


1.38 


1.51 


1.44 


-4.41 


4.2 


4.5 


8 


WN7b 


7.94 


0.71 


0.87 


0.79 


-4.77 


2.4 


6.2 


9 


WN7o 


7.76 


0.78 


0.80 


0.79 


-5.92 


3.8 


5.0 


10 


WN7b 


10.61 


0.65 


0.82 


0.74 


-4.77 


8.5 


3.0 


11 


WN6o 


8.96 


0.73 


0.79 


0.76 


-4.41 


3.3 


5.4 


12 


WN7o 


10.30 


0.82 


0.89 


0.86 


-5.92 


11.8 


6.0 


13 


WC9 


9.45 


0.84 


0.91 


0.87 


-6.30 


9.4 


4.4 


11 


WN6o 


8.71 


1.09 


1.13 


1.11 


-4.41 


2.5 


6.9 


15 


WC8 


9.07 


0.29 


0.47 


0.38 


-4.65 


4.7 


7.0 



crease from 0.55 to 0.90 mag. This increase in extinction 
would reduce our estimated distance from 4.7 kpc to 4.0 kpc. 

Intrinsic absolute Ks-band magnitudes are also seen to 
display a large scatter, even when considering individual WR 
subtypes. With the observed spread in Mj^ s being of the 
order of ±0.5 mag (|Crowther et alj|2006t ). and assuming that 
this dominates the uncertainty in our derived distances, we 
estimate that our heliocentric distances are correct to within 
-25%. 



5 GALACTIC WR DISTRIBUTION 

In order to investigate how the distribution of newly discov- 
ered WR stars compares with that of the previously known 
Galactic WR stars it is necessary to recompute distances to 
the WR stars li sted in the VHth Ca t alogu e of WR stars, 
and its Annex (Ivan der Huchtj l200ll . [2006) using the in- 
trinsic near-IR colour/magnitude calibrations described in 



Section 14.21 For those WR stars which belong to stellar 
clusters or OB associations, we have adop ted previously 
derive d distances from the literature (see Ivan der Huchtj 
(|200lT l; [Crowther et af] |2006) and references therein). The 
only exception is for the highly obscured cluster SGR 1806- 
20, where our derived distance (based on a WN6 calibration) 
suggests a distance of —6 kpc rath er than that of 15 kpc pro- 
posed bv lEikenberrv et al. r <|2004| ). 

For binary members, extinction estimates have been 
made assuming that the WR star dominates the near-IR 
colours of the system. The absolute magnitude of the system 
h as been calculated usin g the synthetic O-star photometry 
of iMartins fc Ple3 (|2006t ). Binarity was only considered for 
cases where the companion star is well classified. 

Estimating the extinction towards dusty WC stars using 
near-IR photometry is not straightforward since the associ- 
ated circumstellar dust can exhibit a wide range of proper- 
ties. At optical wavelengths the contribution of dust is neg- 
ligible so that optically derived extinctions should, in princi- 
ple, be more realistic than those estimated from IR colours. 
Therefore, solely for dusty WC stars, Ak has been estimated 
from optically derived extinctions listed in the VHth Cat - 
alogue assuming that Ay=8.8A/c (|lndebetouw et al1l2005l V 
Distances were then d erived adopting Mk s =-8.5±1.5 mag 
|Crowther et all [20061 '). Applying this calibration to WR95 
yields a heliocentric distance of 3.9 kpc. With the large un- 
certainty associated with Mk s , this consistent with that of 
2.5 kpc inferred from its membership of Trumpler 27. 

In Figure [5] we show the updated Galactic distribution 
of WR stars, as pr ojected onto the Galac tic Plane. We con- 
firm conclusions of (Conti fc Vaccal (Il990h . namely that the 
WR stars appear to trace the spiral structure of the Galaxy, 
with one spiral feature extending towards £~280° . Addition- 
ally, the Galactic WR distribution shows re markable simi- 
larities to that of radio selected Hn regions (|Russeilll2003l . 
their Fig 3). 

The majority of the new WR stars are located between 
4 kpc < Rg < ^©> t ne region where most of the previously 
known WR stars are observed to cluster. The exceptions are 
HDM 10, 12 and 13 which appear to be isolated from the 
bulk of the Galactic WR stars. This distribution suggests 
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Figure 6. The Galactic WN (triangles), WC (squares) and 
WO (circles ) distribution as projected onto the Galactic Plane. 
Open symbols correspond to previously catalogued WR stars 
(van der Hucht|[200ll,l2006h whereas filled symbols represent the 
new WR stars presented here. The Sun (0) corresponds to a 
Galactrocentric distance of 8.5 kpc and the symbol at the centre 
of the plot represents the Galactic Centre, Arches and Quintu- 
plet clusters (N TO t [WR]=60). The 24 WR stars located within 
Westerlund 1 are also represented by a star and have an assumed 
distance of 4.5 kpc. The dashed lines indicates the latitude cover- 
age of the GLIMPSE survey. 



that our survey has identified WR stars situated in more ob- 
scured regions of the inner Galaxy, rather than those which 
are more distant from the Sun. 

Recalling from Section [3, 2 1 arbitrary magnitudes limits 
of Ks —10.5 and Ks — 11.0 mag were placed on IR spectro- 
scopic observations to allow the maximum number of targets 
to be observed. This observational limit obviously influences 
the distribution of the newly discovered WR stars shown in 
Figure |S] and the maximum distance out to which we can 
detect WR stars is, of course, subtype dependant. For ex- 
ample, an early-type WR star (Mk s — -4 mag) at a distance 
of 15 kpc would have a minimum apparent Ks magnitude of 
12 mag, well below the limiting magnitude of any of our fol- 
lowup observations and the GLIMPSE [8.0/im] completeness 
limit of 9.5 mag. 

Adopting an average K-band extinction of 0.86 mag 
(Table[5]) and a limiting magnitude of Ks =10.5 mag, our ob- 
servations are capable of identifying early and late-WN stars 
within —6.5 and 13.5 kpc of the Sun, respectively. Similarly 
for WC stars, observations of WC5-8 and WC9 subtypes are 
limited to within 7.5 and 12 kpc of the Sun, respectively. 

If we were to extend our observations to cover the fainter 
targets in our sample, then it is possible that we could probe 
the WR population at larger distances. However, this may 
prove to be difficult as a combination of larger photometric 
uncertainties and crowding issues may significantly affect 
the reliability of our selection criteria, severely reducing the 
success rate of the survey. 

It is well established that the absolute number and 
subtype distrib ution of a W R population is dependant 
on metallicity (Massev 2003). In the Milky Way, a slight 
metallicity gradient (A log (O/H) = -0.044 ±0.01dex kpc -1 
[Esteban et ail l|2005l )]) is reflected by a preponderance of 



Table 6. Galactic distribution of WR stars as a function of Galac- 
tocentric distance. log(0/H)+12 is quo ted for Rq = 6.0, 8.5 and 
10 kpc, assuming log(Q /H) a +12=8.66 jAsplund et alj|2004l) and 
a metallcity gradient of lEsteban et ail 1 12*0*051) ! 



R G 


<7.0 


7.0-10.0 


>10 


log(0/H) + 12 


-8.8 


8.66 


-8.5 




96 


68 


13 




81 


41 


6 


Nvko 


1 


2 







0.84 


0.60 


0.46 




0.4 


3 


2 


Nwce/Nwcl 


0.01 


1 


2 



late-type WC and early-type WN stars at smaller galactro- 
centric distances. 

In Table [5] we list the number and subtype distribu- 
tion of Galactic WR stars according to our near-IR de- 
rived distances. This shows that the number of WC stars 
relative to WN stars decreases with Galactocentric r adius , 
conf irming the earl ie r con clusions of IConti fc Vaccal (Il990l ) 
and Ivan der HuchtJ l|200ll) . Using oxygen abundance as a 
proxy for metallicity, the early-to-late subtype distribution 
is seen to be extremely sensitive to metallicity, especially 
for WC stars. At Solar metallicity we count an equal num- 
ber of WCE stars relative to WCL stars whereas at higher 
metallicity WCL stars significantly outnumber their early- 
type counterparts. At larger Galactic radii (lower metallic- 
ities) early-type stars dominate the WC population with a 
WCE/WCL ratio of —2. However, with only small number 
statistics, completeness is one important issue that must be 
considered when studying the WR distribution in the outer 
regions of the Milky Way. Similar conclusions are drawn 
for WN stars with late-type stars being more common at 
smaller galactic radii, however the range of observed ratios 
is much less extreme than that for WC stars. These results 
confirm prev ious conclusions drawn about the Galactic WR 
distribution l|Conti fc Vaccalll990l ; Ivan der Huchtll200ll ) but 
highlight how sensitive WR populations are to the metal 
content of their environment. 

Over the last two decades WR populations have 
been w ell samp l ed fo r a variety of metallicty environ- 
ments (iMassevI 120031 ). N(WC)/N(WN) is observed to 
range from —0.1 in the metal-poor SMC (Z~0.2Zq), 
whilst N(WC)/N(WN)-1 for M83. For Solar metal- 
licities, observations of the inner regions of M33 
suggest N(WC)/N(WN)-0.7, very similar to that of 
N(WC)/N(WN)=0.6 observed for 7.0kpc< R G <10kpc. 
For the metal-rich inner Milky Way, Z is equivalent to that 
of M83. In M83, N(WC)/N(WN)-1 comparable to 0.84 ob- 
tained for RG<7.0kpc. 



6 NON-WR DETECTIONS 

In order to better refine our selection criteria, it is vital 
to gain an understanding of the objects that mimic the IR 
colours of WR stars. In this section we will briefly com- 
ment on the objects which satisfied our selection criteria, 
but turned out not to be WR stars. 
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From spectroscopic follow-up observations, ~75% 
of spectra display ed emission features (c.f., 15% in 
iHomeier et al.l 12003), In addition to the 15 WR stars, we 
have identified one previously unreported planetary nebula 
and ~120 (65% of the total sample) emission-line objects. 
Of the 50 spectra (15% of the total sample) which did not 
contain emission features, the majority displayed strong H 
absorption features consistent with mid- to late-type (A-G) 
dwarf stars. 

For the majority (~80 stars) of our emission-line ob- 
jects, emission features consisted of leading H transition 
lines; Ha and H/3 in our optical sample and P/3 and Br7 
in our IR sample. 

In Figure [7J we show the rectified spectra of a rep- 
resentative sample of these objects observed with the 
CTIO/RCSpec. Ha is always present in emission, whereas 
H/3 is seen in emission or absorption. In addition to Ha 
and H/3, the majority of objects display HelA6678 absorp- 
tion/emission. Helium and hydrogen features are charac- 
teristic of early-type stars, but the absence of HellA5411 
is consistent with these objects being early B super- 
giants/hypergiants. Sources which do not show He I A6678 
must be later-type stars, possibly late B or early A super- 
giants. 

We have also identified a subset (~40) of objects which 
display rich emission-line spectra, including H, He I and nu- 
merous lower-excitation metal lines, especially Fell. Objects 
which display these emission-line features and near-IR ex- 
cesses are usually classified as peculiar B[e] stars. A compari- 
son of one these objects with that of the B[e] star HD 326823 
is presented in Fig [7] clearly illustrating the similar spectral 
morphologies. Consequently, based solely from a spectral 
comparison standpoint, we identify these objects as Be/B[e] 
candidates. 



7 DISCUSSION & CONCLUSIONS 

Using a combination of near- and mid-IR colours we have 
identified candidate WR stars in the Galactic Plane using 
the observed colours of known Galactic WR stars as a "train- 
ing set" . Initial spectroscopic follow-up of 184 of our candi- 
date WR stars has led to the discovery of eleven new WN 
and four new WC stars. 

Up to now, WR searches have used narrowband K-band 
interference filters to search for emission line excesses, mak- 
ing them particularly insensitive to dusty WC stars. Here, 
we have detected both WN and WC stars, with none of the 
WC stars showing evidence for circumstellar dust. 

The fifteen WR stars fall within 10°. 2 < I < 331°. 0. 
For fourteen of these, we estimate that Rg < R© suggesting 
that this survey has been sampling WR stars located in 
nearby, highly obscured regions of the Galaxy. Since an all 
sky survey was not possible and follow-up spectroscopy was 
only capable of detecting the faintest WR subtypes to within 
~6.5 kpc, this survey is helping to complete the WR census 
within the Solar vicinity. 

Many recent WR surveys have concentrated on find- 
ing WR stars in massive stellar clusters and with 30% 
of th e Galactic WR stars being contained with four clus- 
ters (|van der HuchtJ 120061 ') . such searches are very worth- 
while. Interestingly however, eight of the nine WR stars de- 
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Figure 7. Non-WR detections. Top: A montage of CTIO/optical 
spectra of the Ha emission-line objects detected. Bottom: A 
spectral comparison of 2MASS source 11585474-6234116 with 
HD 326823 (B[e]). The spectra have been been vertically offset 
(by 1.5 continuum units) for clarity. Prominent emission features 
are marked. 



tected in recent pho tometric surveys (jHopewell et alj|2005l : 
IHomeier et al.ll2003l ) and the fifteen WR stars presented here 
are located in relatively isolated regions. Continually find- 
ing WR stars in isolation raises many questions as to how 
they arrive there, are they runaway stars expelled from clus- 
ters by supernovae events or binary/close interactions, have 
their native clusters simply dispersed or do some massive 
stars actually form in isolation? Of course, with only small 
number statistics we cannot attempt to address such issues 
here, but it does highlight the importance of sample com- 
pleteness and that searches for WR stars in the Galactic 
Plane are extremely worthwhile. 

In addition to WR stars, we have identified one new 
planetary nebula and ~120 H emitting sources within the 
Galaxy, most of which we believe to be massive OB su- 
pergiants/hypergiants. With a WR yield rate of ~7%, and 
an massive star detection rate of ~65%, our initial findings 
have proved that this survey is one of the most successful 
searches for evolved massive stars in the Galaxy. However, 
since GLIMPSE is restricted to >±10° from the Galctic 
center, it is not possible that we will find all of the missing 
WR stars, but will, in fact, complement other ongoing WR 
searches. 
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Our survey has revealed that using a combination of 
2MASS and GLIMPSE data it is possible to locate WR 
stars which had remained hidden from view behind copi- 
ous amounts of dust. Assuming that subsequent follow-up 
observations of our remaining candidates achieves the same 
success rate, we have the possibility of increasing the known 
Galactic WR population twofold. 
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Figure Al. Optical CTIO/RCSpec data for six of the newly discovered WR stars. WN stars are given in the top panel and WC star 
are presented in the bottom panel. 
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Figure A2. Ncar-IR spectra of five of the newly confirmed Galactic WR stars. WC stars are shown in the top panel and WN stars are 
bottom panel. The spectra were acquired using IRTF/SpeX, except for HDM 14 which was observed with APO/CorMASS. 
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Figure A3. Near-IR spectra of the four newly confirmed WN7 stars. Spectra have been grouped into broad (top panel) and narrow-lined 
(bottom panel) stars. The spectra were acquired using IRTF/SpeX, except for HDM 12 which was observed with APO/CorMASS. 



14 Lucy. J. Hadfield et al. 




HDM 2 **, j. 

9' . 
II 
1 

* • 
• > 

• 4 


HDM 3 * » 

, | 

• ■ t m 

i * i 

; / ■ 


1 s ■ . 

HDM 5 

* B 

1 ■ 

• ♦ 

m , % v * 
l* * 

f A 


- hdm ^ * • y 

* Jf ' * ■ 

. m ... 


HDM 8, ' • \ 

v. 

. • » « 

M * • 


Hrtj?9*;" ■ * * 

• * • 


HDMllV - 

• v 

• 1 . » • • 
• 

«r * 


»" • V* 

■HDM 14 • » • 
■ ■ 

•; H i ^ 


HDM 14 • „ 
■ Ml * • 

• : -4^ * ' 

► " • «• • • • * . . . 

• • 

• " * \ 


' HDM 15, 

* .0* 

.1 

* • % 4 



Figure Bl. 2x2' 2MASS Kg finding charts for the newly discovered WR stars. North is up and east is to the left on all images. 



